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[摘要 ] 2019 年 12 月 ， 新 型 冠状 病毒 肺炎 (NCP， 又 称 2019-nCoV) 疫 情 从 武汉 开始 爆发 , 几 天 
内 迅速 传播 到 全 国 乃至 海外 ， 对 我 国 的 工农 业 生 产 和 人 民生 活 产生 了 重要 影响 。 科 学 有 效 
掌控 疫情 发 展 对 疫情 防 控 至 关 重 要 。 本 文 基于 中 国 卫 健 委 及 湖北 省 卫 健 委 每 日 公布 的 累计 
依 数 ， 采 用 逻辑 斯 带 模 型 对 数据 进行 了 拟 合 ， 以 期 给 该 疾病 的 防 控 治 提供 科学 依据 。 通 

过 公布 的 疫情 数据 ， 我 们 反 演 了 模型 的 参数 ， 进而 有 效 地 模拟 了 目前 疫情 的 发 展 ， 并 预测 
了 疫情 未 来 的 趋势 、 我 们 预测 ， 湖 北 省 疫情 还 要 持续 至 少 2 周 ， 而 在 全 国 其 他 地 区 ， 疫 情 
可 望 1 周 左右 达 到 顶峰 。 
[关键 词 ] 冠状 病毒 感染 肺炎 数学 预测 湖北 省 全 国 其 它 地 区 
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Abstract 


Novel Coronavirus Pneumonia (NCP, or alternatively 2019-nCoV), initially blown up 
in Wuhan in December of 2019, has been quickly spread all over China, and even 
other countries of the world, which has produced an Important effect to the 
agricultural and industrial activities, and daily life. It is expected that a well-known 
recognition ls essential to the effective prevention of the disease. Based on the daily 
announced numbers ofthe infective 


people from the National and Hubei provincial Health commissions, a logistic model 
is applied in this paper for data fitting, in order to provide some scientific information 
for the effective prevention and controlling of the disease. Using the parameters 
obtained from the data simulation, a susceptible-infected (SI) model is used to 
forecast the future trend of the NCP. Our work indicated that the epidemic will last at 
least two additional weeks in Hubei, but should come to an apex in one week in other 
areas of China. 
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1 引言 

2019 年 底 至 2020 年 初 ， 我 国武 汉 爆 发 新 型 冠状 病毒 感染 的 肺炎 ， 引 起 全 国人 民 的 关 
自 2020 年 1 月 23 日 武汉 封 城 至 今 已 经 三 周 ， 该 传染 性 疾病 的 发 展 趋势 对 未 来 如 何 部 
所 新 的 扩 近 泊 渍 放 乃至 国内 生产 的 恢复 都 有 着 非常 重要 的 影响 。 鉴 于 这 些 与 人 民生 活 密 


切 相关 的 重要 问题 ， 国 内 外 学 者 运用 不 同 的 数学 模型 ， 对 该 疾病 的 扩散 与 传播 做 了 一 定 的 
工作 [1-3]， 但 就 该 疾病 未 来 的 发 展 趋势 ， 目 前 开展 的 研究 并 不 多 见 。 因 此 ， 本 文 利 用 目前 
被 大 多 数学 者 认可 的 数学 模型 ， 对 冠状 病毒 感染 肺炎 疾病 的 发 展 趋势 进行 预测 ， 以 期 对 目 
前 的 防 控 治 工作 提供 科学 参考 。 


2 模型 建立 

设 y 为 已 感染 人 数 ，t 为 时 间 ， 在 总 人 数 趋 于 无 穷 的 情况 下 (或 者 说 ， 在 已 感染 人 数 
相对 总 人 数 是 一 个 可 以 忽略 的 量 的 情况 下 ) ， 此 时 ， 感 染 人 数 增加 的 速率 与 已 感染 人 数 成 
正比 ， 因 此 我 们 可 以 得 到 常 微分 方程 y”=ay。 解 这 个 微分 方程 ， 我 们 可 以 得 到 ， 等 比 数列 
上 升 的 趋势 图 ， 这 也 是 疫情 初期 感染 人 数 差不多 是 等 比 上 升 的 原因 。 

不 过 ， 总 人 数 其 实 是 有 限 的 。 尤 其 是 在 做 了 隔离 措施 后 ， 一 个 人 能 接触 的 人 也 是 有 限 
的 。 这 样 ， 感 染 人 数 的 增加 实际 上 与 已 感染 人 数 和 未 感染 人 数 乘 积 成 正比 。 再 加 上 有 了 隔 
离 措施 ， 切 断 了 大 多 数 人 与 病人 的 联系 ， 因 此 我 们 可 以 列 出 如 下 常 微分 方程 y”=ay(b- 
y) 。 其 中 ，a 代表 传染 率 ，b 代表 可 能 感染 的 总 人 数 。 这 个 方程 的 解 要 稍微 复杂 一 些 ， 可 以 
通过 以 下 步骤 解 出 : 


那么 最 终 ， 我 们 可 以 得 到 

实际 上 ， 这 个 模型 就 是 传染 病 的 SI 模型 [4] 。 其 微分 方程 的 解 有 解析 表达 式 ， 易 于 优 
化 求解 ， 但 没有 考虑 病人 的 康复 ， 以 及 潜伏 期 的 情况 。 若 考虑 康复 ， 则 得 到 SIR 模型 
[5] ， 若 再 考虑 潜伏 期 ， 即 得 到 SEIR 模型 ， 甚 至 更 复杂 的 SEIJR 模型 [6] 。 尽 管 传染 病 研究 
多 采用 SIR 或 SEIR 模型 ， 但 这 两 个 模型 的 解 均 没有 显 式 表达 ， 这 意味 着 需要 用 数值 方法 
解 常 微分 方程 ， 这 意味 着 后 面 的 数据 拟 合 ， 也 就 是 优化 过 程 将 会 非常 复杂 ， 且 在 疫情 初期 
康复 人 数 相 比 感染 人 数 可 以 忽略 不 计 。 因 此 ， 我 们 仍 采用 相对 简单 的 SI 模型 ， 以 期 得 到 较 
为 可 靠 的 模拟 结果 。 
这 些 参数 中 ，A 代表 传染 的 速率 ，k 代表 最 终 传染 病人 累计 的 数量 ，K/P 大 概 是 第 一 天 
新 型 冠状 病毒 肺炎 的 病例 数量 。 


3 数据 拟 合 

上 述 讨论 只 是 一 个 理论 分 析 ， 并 没有 结合 本 次 疫情 的 数据 。 而 且 ， 其 中 的 k, b, C, a 等 
变量 都 极 难 通 过 经 验 确 定 ， 因 此 我 们 要 结合 我 们 已 知 的 数据 ， 确 定 那 些 变 量 以 及 K,P,A 的 
值 ， 这 个 过 程 就 是 数据 拟 合 的 过 程 。 


1 冠状 病毒 肺炎 疫情 


表 
增 诊 增 
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31161 |3143 |22112 |2447 |9049 |696 


全 国情 况 : 这 里 ，N(1)=214 是 截至 2020 年 1 月 19 日 24 时 ， 全 国 累计 报告 新 型 冠状 病 
毒 感染 的 肺炎 病例 数量 。N(20)= 34546 是 截至 2020 年 2 月 7 日 24 时 ， 全 国 累计 报告 新 型 
冠状 病毒 感染 的 肺炎 病例 数量 。 我 们 尝试 根据 这 20 个 点 拟 合 出 曲线 ， 并 根据 拟 合 出 的 曲线 ， 
对 接 下 来 7 天 的 全 国 确诊 病例 数量 进行 预测 。 这 样 ， 我 们 的 任务 就 是 求解 如 下 最 优化 问题 


此 优化 问题 几 个 变量 尺度 差距 非常 大 ， 是 一 个 高 度 病态 问题 ， 好 在 维 数 不 高 ， 可 以 通 
过 和 牛顿 法 或 高 斯 牛顿 法 求解 。 我 用 高 斯 牛顿 法 求 得 下 图 (图 1) 。 
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图 1 全 国 感染 新 型 冠状 病毒 肺炎 情况 


在 该 图 中 ， 圈 是 用 来 构建 模型 的 训练 数据 ， 又 是 训练 数据 后 7 天 的 真实 结果 ， 用 于 验 
证 预测 的 可 靠 性 。 从 图 中 可 以 看 出 ， 这 条 曲线 的 预测 并 不 可 靠 。 它 虽然 很 好 地 拟 合 了 2 月 
日 之 前 的 状况 ， 但 后 面 的 预测 误差 较 大 。 

此 外 在 发 病 时 期 ， 湖 北 省 与 国内 其 它 地 区 在 医疗 条 件 等 方面 有 明显 不 同 ， 因 此 本 文 将 
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数据 ， 但 给 出 的 预测 与 实际 有 较 大 差距 。 
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图 3 全 国 其 它 地 区 感染 新 型 冠状 病毒 肺炎 情况 


2 月 17 日 2 月 22 日 2 月 27 日 


根据 上 述 同样 的 原理 ， 我 们 得 到 图 2 和 图 3。 尽管 它们 都 很 好 地 拟 合 了 2 月 7 日 之 前 的 


我 们 猜想 ， 这 是 由 于 1 月 27 日 ， 新 的 测试 试剂 开始 应 月 
从 而 导致 确诊 人 数 大 幅 提升 。 或 者 说 ， 可 以 


较 大 的 差距 。 因 此 ， 我 们 又 以 1 月 27 日 


说 ，26 日 及 以 前 


日 ， 而 这 
『， 所 


以 来 的 数据 重 刘 


文 导致 确诊 率 的 大 幅 上 升 ， 
的 数据 都 和 真实 感染 人 


进行 模拟 ， 结 果 如 图 4 图 5 所 
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4 湖北 省 感染 冠状 病毒 肺炎 情况 (1 月 27 日 始 ) 


该 组 数据 不 仅 很 好 地 拟 合 了 2 月 7 日 之 前 的 数据 ， 还 对 后 7 天 的 数据 提供 了 很 好 的 佑 
计 ， 拟 合 效果 较 好 。 拟 合 得 到 最 终 湖北 省 感染 新 型 冠状 病毒 肺炎 人 数 累 计 应 在 40730 人 左 
右 。 我 们 认为 ， 这 可 能 是 一 个 比较 保守 的 估计 ， 真 实 人 数 要 密切 根据 动态 变化 情况 予以 调 
整 。 其 实 ，2 月 5 日 以 后 ， 每 天 湖北 省 确诊 人 数 的 增加 量 都 较 前 日 有 所 降低 ， 可 以 确认 拐 
点 已 经 过 去 了 。 

再 考虑 全 国 其 它 地 区 的 情况 ， 过 程 与 上 图 类 似 。 模 拟 估计 最 终 全 国 除 湖北 省 其 他 地 区 
感染 新 型 冠状 病毒 肺炎 人 数 约 为 11630 人 ， 但 现在 感染 人 数 已 经 超过 了 这 个 数字 。 这 条 曲 
线 对 后 面 的 几 天 拟 合 的 依旧 不 好 ， 对 此 ， 本 文 认为 可 能 是 由 以 下 几 点 原因 导致 的 : 

1， 模型 不 够 精确 。 

本 文采 用 了 SI 模型 ， 这 个 模型 在 疫情 初期 是 对 复杂 的 SEIR 模型 很 好 的 近似 。 但 是 现 
在 疫情 已 经 跨 过 了 拐点 很 多 ， 并 且 治 愈 人 数 已 经 有 了 8000 多 人 ， 这 对 于 总 确诊 病例 近 5 万 
人 来 说 已 经 是 不 可 忽略 的 量 了 。 

2， 确诊 率 在 提高 。 

随 着 时 间 的 推移 ， 可 能 医学 研究 的 发 展 ， 医 疗 物资 的 状况 都 有 了 很 大 的 变化 ， 确 诊 率 
一 定 会 随 着 时 间 不 断 提 高 。 

3. 数据 不 够 准确 ， 有 误差 。 

统计 数据 时 由 于 重复 ， 漏 报 等 因素 ， 会 导致 确诊 病例 与 真实 感染 人 数 有 一 定 差 据 ， 这 
些 差距 在 数据 拟 合 的 过 程 中 会 放大 。 
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图 5 全 国 其 它 地 区 感染 冠状 病毒 肺炎 情况 (1 月 27 日 始 ) 


用 1 月 27 日 到 2 月 7 日 的 全 国 感染 数据 ， 我 们 可 以 得 到 图 6. 这 个 预测 还 是 比较 可 靠 的 ， 从 
图 中 ， 我 们 可 以 看 到 ， 最 终 全 国 感染 冠状 病毒 肺炎 累计 人 数 大 约 在 5.3 万 人 左右 。 
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图 6 全 国 感染 冠状 病毒 肺炎 情况 (1 月 27 日 始 ) 


总 之 ， 根 据 拟 合 与 优化 模型 计算 ， 湖 北 省 感染 新 型 冠状 病毒 肺炎 的 人 数 累 计 应 在 4. 1 
万 人 左右 ， 全 国 其 它 地 区 在 1.3 万 人 左右 ， 共 计 约 5.4 万 人 。 未 来 ， 可 根据 数据 实时 更 新 
拟 合 结果 ， 以 提供 较 可 信和 的 数据 。 


下 面 我 们 讨论 本 次 疫情 何 时 达到 拐点 的 问题 。 所 谓 拐点 就 是 新 增 病例 数 不 再 增加 。 表 
1 的 实际 数据 显示 ， 湖 北 省 新 增 病例 在 2 月 4 日 达到 高 峰 ， 但 随后 数 日 一 直 处 于 波动 状态 ， 
直至 2 月 9 日 以 后 ， 才 基本 呈 下 降 趋 势 。 与 湖北 省 的 情况 不 同 ， 全 国 其 它 地 区 自 2 月 3 日 
以 来 ， 新 增 病例 已 经 连续 下 降 ， 应 该 标明 拐点 早已 经 到 来 。 
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图 7 全 国 新 增 感染 


尽管 我 们 的 模拟 计算 表明 疫情 拐点 已 过 (图 7)〉 ， 但 我 们 仍 需 提 


措施 不 到 位 ， 


导致 更 多 的 人 被 传染 ， 从 而 爆发 第 二 波 疫情 。 


本 次 疫情 何 时 结束 也 是 大 家 比较 关心 的 话题 。 0 
增 病例 降 至 10 为 参考 点 ， 湖 北 省 将 持续 至 3 月 1 日 ， 全 国 其 它 


就 是 说 ， 如 果 


拉 冠 状 病毒 肺炎 拟 合 情 况 (1 月 27 日 始 ) 


高 警惕 ， 


地 区 将 持续 


3 月 1 日 


至 2 月 25 日 。 


保持 目前 防 控 措 施 的 话 ， 本 次 疫情 大 约 在 3 月 初 后 将 得 到 根本 控制 。 


本 文采 用 
毒 肺炎 人 数 大 
值得 指出 


目前 广 为 通 用 的 SI 传染 病 模型 ， 


过 拟 合 和 优化 计算 ， 获 得 全 国 感染 


我 们 以 新 


3 月 6 日 


防止 因 防 控 


也 


冠状 病 


er eta Ve tes Bt 全 国 其 它 地 区 1. 3 万 人 左右 。 


的 ， 我 们 的 模型 对 全 国 其 他 地 区 的 情况 拟 合 得 较 差 ， 可 能 需要 采 上 月 


SEIR 模型 ， 但 该 模型 参数 求解 难度 较 大 。 男 一 方面 ， 就 本 文采 | 


加 大 近期 权重 


2 月 13 日 ， 


口径 的 变动 对 


天 来 ， 由 此 而 带 来 的 湖北 省 
大 于 我 们 预测 的 总 感 


本 文 预 测 数据 
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的 模型 而 言 ， 也 可 以 尝 
这 些 将 在 下 一 步 工 作 中 尝试 。 
新 冠 肺炎 确诊 也 可 以 包括 临床 确诊 方法 。 
但 将 显著 提升 湖北 省 感染 人 数 。3 
肺炎 感染 人 数 上 升 17400 余 人， 感染 总 人 数 达 56000 余 人 , 已 
染 人 数 。 如 果 这 一 数据 准确 可 靠 的 话 ， 湖 北 省 的 实际 感染 人 数 将 超出 


尺 乱入 
AN [2 多 一 
~ 请 ) 


40% 左 右 ， 且 疫情 结束 的 日 期 将 推 后 两 周 左右 ， 值 得 相关 部 门 重视 。 
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导 下 完成 ， 


导师 责 亚 湘 研究 员 指 并 得 到 国家 
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